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ABSTRACT 

The source of Faraday rotation in the jet of the radio galaxy 3C 120 is analyzed through Very 
Long Baseline Array observations carried out between 1999 and 2007 at 86, 43, 22, 15, 12, 8, 5, 2, 
and 1.7 GHz. Comparison of observations from 1999 to 2001 reveals uncorrelated changes in the 
linear polarization of the underlying jet emission and the Faraday rotation screen: while the rotation 
measure (RM) remains constant between approximately 2 and 5 mas from the core, the RM- corrected 
electric vector position angles (EVPAs) of two superluminal components are rotated by almost 90° 
when compared to other components moving through similar jet locations. On the other hand, the 
innermost 2 mas experiences a significant change in RM - including a sign reversal - but without 
variations in the RM-corrected EVPAs. Similarly, observations in 2007 reveal a double sign reversal 
in RM along the jet, while the RM-corrected EVPAs remain perpendicular to the jet axis. Although 
the observed coherent structure and gradient of the RM along the jet supports the idea that the 
Faraday rotation is produced by a sheath of thermal electrons that surrounds the emitting jet, the 
uncorrelated changes in the RM and RM-corrected EVPAs indicate that the emitting jet and the 
source of Faraday rotation are not closely connected physically and have different configurations for 
the magnetic field and/or kinematical properties. Furthermore, the existence of a region of enhanced 
RM whose properties remain constant over three years requires a localized source of Faraday rotation, 
favoring a model in which a significant fraction of the RM originates in foreground clouds. 
Subject headings: galaxies: active - galaxies: individual (3C 120) - galaxies: jets - polarization - radio 
continuum: galaxies 



1. INTRODUCTION 

Polarimetric Very Long Baseline Interferometric 
(VLBI) observations have proven to be a very power- 
ful tool to probe the magnetic field structure in rela- 
tivistic jets of Active Galactic Nuclei (AGN). However, 
recent observations have revealed that the actual orien- 
tation of the polarization angle at parsec scales m ay be 
severely affected by Faraday rotation (e.g., Udomprasert] 
eFaL] [19971 |Tay lor | [19 98; Reynolds et al. 2001; Asada 



et al."20 02[|Ziavala 
2005; Gomez et al. 



Kharb et al.||2009 



^ Taylor,2003, .2004, .Attrid ge et al. 
20081 lo'Sullivan Gabuzd a 2009; 



Asada et al.||2010| [Taylor Zavala 



2010), which rotates the polarization vector by an angle 
proportional to the square of the observing wavelength. 
The analysis of Faraday rotation in jets of AGN is there- 
fore of special relevance to properly determine the actual 
orientation of the polarization angle, and hence of the 
magnetic field structure on the plane of the sky in the 
emitting region. 

Several possibilities have been suggested as the source 
of the observed Faraday rotation measure (RM) in jets 



of AGN. Gradients in the RM along the jet (e.^., [Taylor 
T9981 jGabuzda fc Chernetskii|[2QQ3l [Gomez et alj^^ 



and variabili ty in both s pace an d time (e.gT^ Zavala <fe 
p:aylor 20^ [20031 [20051 |Asada eFari[2Q08a| |Mahmu5 
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et al.|[2QQ9 ) suggest that the observed RMs are intrinsic 
to the mner few hundred par sees of AGN. This rules 
out the possibility that the RM is produced by our own 
Galaxy, or in the host galaxy interstellar or intracluster 
medium. Internal Faraday rotation i n the emitt i ng jet 
would produce severe depolarization ( Burn 1966| [Cioffi 



& Jones 1980) that is not observed, and can therefore 
be safel y excluded as the cause of RM in most c ases (see 



Gomez et alT1[2QQ8l [Homan et aLl[2QQ9l ) 



however 

The broad emission-line region (BLK) is thought to be 
smaller than a parsec and to have a small volume cover- 
ing factor, so it can therefore be eliminated as the source 
of the RM that appears on scales of tens or hundreds of 
parsecs. The narrow emission-line region (NLR) extends 
to hundreds of parsecs, and could therefore in principle 
contribute to the observe RM, although its relevance is 
subject to uncertaint ies owing to its largel y unknown vol- 
ume covering factor ( Netzer fc Laor|1993 ). The observed 



RM time variability has also been used as an argument 
to exclude the NL R, or the intercloud gas, as the main 
source for the RM (Asada et al |2008a| [Kharb et al. 12009]) . 

First suggestion that the KM is produced in a sheath 
that is in clo se proximity with the jet was given by [Asada 



et al. (2002), based on the detection of a RM gradient 
across the jet width in 3C 273. Such transverse gradi- 
ents are expected if the sheath is threaded by a helical 
magnetic field, since they can be produced by a system- 
atic change in the net line-of-sight magnetic field com- 
ponent across the jet (Laing 1981; Blandford 1993). Ob- 
servational confirmation for the existence of these helical 
magnetic fields is of special relevance, since it is thought 
that these may appear naturally through the rotation of 
the accretion disk or black-hole ergosphere from which 



2 



Gomez et al. 



jets are launched, and could have an important role in 
the actual f ormation, collimat i on, and acceleration pro- 
cesses (e.g.,[K oide et aL"2QQ2'; 'Vlahakis fc Konigl"2QQ4! 
de Vihiers et al. 2005; Hawley &; Krolik 2006; McKinney 
>06 IKo rmssarov et al. 2007; McKinney & Blandfordj 
>Q9l [Marscher et al. 2008 , 2010; Jors tad et al. 2010^ 
Further observatioins have conhrmed the transverse 
RM gradient in 3C 273, revealing also time variability 
on a s cale of years (IZavala & Taylor||2005| lAsada et al. 



2008a). Following this initial detection m 3 (J 273, other 



sources have also been claimed to contain transverse gra- 
dients h]__th^e_^M al. 2004; Gomez et all 
20081 lAsada et al.n200^bi lO^^ullivan Gabuzda 2009; 



Kharb et al.][ 2Q09; Croke et al. 2009| [AsadaTT al. 2010). 



Contopoulos et al.| ( |200 9) present a study of the observed 
KM in jets of AGNT^laiming that 29 out of 36 sources 
possess transverse gradients in RM, further suggesting 
that the gradi ents have a preferred s ense of direction. 
More recently, Taylor & Zavala (2010) have raised con- 
cerns about the claimed transverse KM gradients by ar- 
guing that (besides other reasoning based on jet opacity 
and significance of the measurements) most of these ob- 
servations lack the necessary resolution transverse to the 
jet. iTaylor fc Zavalaj ( |2010[ ) propose several criteria for 
establishing a reliable transverse RM gradient, claiming 
that only observations of 3C 273 have met these criteria 
so far. 

The existence of a sheath that surrounds a fast, 
synchrotron-emitting j et is also s upported by multiple 
numerical simulations (lAloy e t al.|1999. 2000; Tsinganos 



Bogovalov 2002; de Villiers jet al 
2005 



racia et al 



McKinney 2006 



,2009^ iHawley fc Krohk| |2006| ^ 

ermore, computation of the expected rotation mea- 
sure from 3D gen eral relativistic m agnetohydrodynamic 
simulations by Broderick fc McKinney (2010) show that 
it is possible to reproduce many of the KM structures 
in jets of AGN, including the presence of transverse RM 
gradients. There is, therefore, observational and theoret- 
ical support for the existence of a sheath that surrounds 
the emitting jet that may contribute significantly to the 
observed RM in jets of AGN. 

Addition ally, there is also st rong observational evi- 
dence (e.g., Oosterloo et al. 200Q{ [Solorzano-Iharrea etal 



2001[ IMiddelberg et al. 2007) an d theoretical support 
" "elg:, Stetfen et al. 1997b a; Wang et al.| |2000| |BickneII| 



et al. [2003( [Saxton et al. 1 ,2005) for the existence of tore- 
ground clouds, probably mter acting or entrained by the 
j et, that may contribute significantly to the o bserved RM 
([Gomez et al. 2008 Mantovani et al. 2010). In partic- 



ular, Gomez et al. ( 2008p present a study of the radio 
galaxy 3G 120 consistmg of 12 monthly polar imetric 15, 
22, and 43 GHz Very Long Baseline Array (VLBA) ob- 
servations revealing the existence of a localized region 
of high (~ 6000 rad m~^) Faraday rotation measure be- 
tween approximately 3 and 4 mas from the core. The 
detection of this region of high RM required the use 
of a novel method, consisting of the combination of the 
RM information gathered at all epochs, so that multiple 
superluminal components can be used to map the RM 
screen. This technique revealed a remarkably consistent 
RM across epochs. A smooth sheath around the jet in 
3C 120 cannot produce a localized region of enhanced 
RM - although it successfully explains the observed gra- 
dients in RM and degree of polarization along and across 



the jet - so the RM was assumed to originate in a fore- 
ground cloud, presumably interacting with the jet. 

In this paper we present further observations of the ra- 
dio galaxy 3C 120 aimed to obtain a better understand- 
ing of the origin of the RM screen in this source. Thanks 
to its proximity (z=0.033), VLBI observations are able to 
resolve the jet across its width, revealing a very rich and 
dynamic structu re, even when obse rved at the highest 



frequencies (e.g., [Gomez et al.||1998l [199 9, 2000, 200Tbl 
Walker et al||1987l|2001[ [Jorstad etar[|2005( iMarscKer 
et al 20021 [2007[ [Chatterjee et al.[[2009| [Roc~ bogorb 



et al.[[2QTQ[ . 3C120 is located at a distance of 140 Mpc 
{Hq = 71 km s~^ Mpc~^), so that 1 mas corresponds 
to 0.64 pc. Superluminal motion of components along 
the jet provide an uppe r lim it of the jet viewing angle of 
~ 20° (Gomez et al. 2000[ ), so that distances along the 
jet appear appreciably foreshortened by projection. 

2. OBSERVATIONS AND IMAGES 
2.1. Observations and data reduction 

We present VLBA observations of the jet in the radio 
galaxy 3C 120 corresponding to two different observing 
programs carried out in 2001 and 2007. Partial anal- 
ysis of these data re vealed gradients in Faraday rota- 
tion and polarization ( Gomez et al.|20Q8 ), and un usually 
high br ightness temperature 140 pc from the core ( Roca- 
[Sogorb et al. 2010). 

Observations during 2001 were made in dual polariza- 
tion at the standard frequencies of 15, 22, and 43 GHz, 
covering a total of 12 epochs: 2000 December 30 (here- 
after epoch A), 2001 February 1 (B), 2001 March 5 (C), 
2001 April 1 (D), 2001 May 7 (E), 2001 June 8 (F), 2001 
july 7 (G), 2001 August 9 (H), 2001 September 9 (I), 
2001 October 11 (J), 2001 November 10 (K), and 2001 
December 13 (L). The data were recorded in 2-bit sam- 
pling VLBA format with 32 MHz bandwidth per circular 
polarization at 256 Mbits/s, except for epochs B, C, D, 
F, and G, for which 1-bit sampling and 128 Mbits/s were 
used. One of the Very Large Array (VLA) antennas was 
added to the array at epochs A and E. Bad weather con- 
ditions and/or hardware problems, resulting in the loss 
of data, occurred at: HN during epochs A, C and D; OV 
for epoch H; and KP, MK and BR at epoch L. 

The 2007 observations were performed on November 
3 (at 86 GHz), November 7 (43, 22, and 15 GHz), and 
November 30 (12, 8, 5, 2, and 1.7 GHz) with ah available 
VLBA antennas except Saint Croix, which was down for 
maintenance. The data were recorded in 2-bit sampling 
VLBA format with 32 MHz bandwidth per circular po- 
larization at a recording rate of 256 Mbits/s. Observa- 
tions at 12 GHz were split into two 16 MHz bandwidths 
centered at 12.59 and 12.12 GHz to maximize possible 
detection of small Faraday rotation on relatively small 
angular scales. 

Reduction of the data was perf ormed with the AIPS 
software in the usual manner (e.g., Leppanen, Zensus, & 
Diamond 1995"). Opacity corrections at 22, 43, and 86 
GHz were introduced by solving for receiver temperature 
and zenith opacity at each antenna. The feed D-terms 
(instrumental polarization) were found to be very consis- 
tent over all observed sources and to remain stable across 
epochs during the 2001 observations. 

The absolute phase offset between the right- and left- 
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Fig. 1. — 15 (left), 22 (middle), and 43 GHz (right) natural weighted VLBA images of 3C 120 in 2001. Note the different scale size used 
for each frequency. Vertical map separations are proportional to the time difference between epochs, shown in 15 GHz panel. Total intensity 
contours are overlaid at 0.5 (0.8; 1.0), 1.1 (1.7; 2.2), 2.5 (3.6; 4.6), 5.6 (7.7; 9.9), 13 (16; 21), 28 (35; 45), 62 (74; 95), 138 (158; 203), 310 
(337; 433), and 692 (716; 923) mJy beam"! at 15 GHz (22; 43 GHz). A common convolving beam of FWHM 1.14 x 0.54 (0.78 x 0.37; 
0.40 X 0.19) mas at -1.3° (-0.9°; -1.1°) at 15 GHz (22; 43 GHz) was used for all images and is shown in the lower left corner of each image 
sequence. Color images (on a linear scale shown at the top of each sequence) show the linearly polarized intensity. White bars (of unit 
length) indicate the electric vector position angle, uncorrected for Faraday rotation. 



circularly polarized data, which determines the electric 
vector position angle (EVPA), was obtained by compari- 
son of the integrated polarization of the VLBA images of 
several calibrators (0420-014, OJ 287, BL Lac, 3C 454.3, 
DA193, 3C 279, and 4C 39.25) with VLA observations at 
epochs 2000 December 31, 2001 February 3, 2001 March 
4, 2001 April 6, 2001 May 11, 2001 June 10, 2001 August 
12, 2001 September 10, 2001 October 12, 2001 November 
6, 2001 December 15, and 2007 November 19, as well as 



archival data from the UMRAO, MOJAVE, and NRAO 
long term monitoring programs. Estimated errors in the 
orientation of the EVPAs vary from epoch to epoch and 
with frequency, but usually lie in the range of 5°-10°. 
After the initial reduction, the data were edited, self- 
calibrated, and imaged both in total and polari zed inten- 
sity w ith a combination of AIPS and DIFMAP ( [Shepherd 
1997|). Comparison of the D-terms across epoc hs pro^ 



vides an alternative calibration of the EVPAs (Gomez 
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Fig. 2.— VL BA image s of 3C 12 in 2007 November 30 (2007.91 ) at 1.7, 2, 

intensitv (naturally weighted) contours are overlaid at 0.9, 2.0, 4.2, 8.9, 19, 41, 86, 

oTi^'WHM 16.85x5.55, 11.85x4.11, 5.8x1.9, and 3.36x1.11 mas at -18^ tor frequencies 1.7, 2, 5, and 8 GMz^^esp ect i vel v.|^^olors_show_ttie 
Fmearly polarized intensity, and bars (of unit length) indicate the EVPAs, uncorrected for Faraday rotation, 



) and 8 GHz , from t op to bottom, r espectively. Total 
184 , 392, and 834 mJ y beam~^, with c onvolving beams 
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Fig. 3. — Same as Fig. [2I for observati ons at 12 (2 IFs centered at 12.59 GHz), 15, 22, and 43 GHz, from top to bottom, respectively. 
Ubserv afions were fallen m 2UU7 INovember Y (:^U UY.85), except at 12 (JHz (epoch 2UUY.91). (Jontours are overlaid at 1.1, 2.4, 5.3, 11, 25, 
53, 115 , 249, 538, and 1164 mJy beam-\ with c onvolving beams of FWHM 2.21x0.71, 1.76x0.61, 1.25x0.4, and 0.63x0.22 mas at -18° 
or frequencies 12, 15, 22, and 43 GHz, respectively. 



et al.||2QQ2 |, which was found to be consistent with that 
obtained by comparison with the VLA data. 

2.2. 2001 images 

Figure [1] shows the total and linearly polarized inten- 
sity images corresponding to the 2001 monitoring pro- 
gram. These were used to obtain the sequence and mean 
value of the rotation measure maps shown in Figs. 1 and 
2 of ^Gomez et al.^ (2008^ , respectively, which requires 



careful alignment of the images across epochs and fre- 
quencies. A detailed description of the procedure used 
and accuracy achieved is presented in Appendix [A) 

2.3. 2007 images 

VLBA observations of 3C 120 taken during November 
2007 are shown in Figs. |2j [3] and [4j which provide a 50- 
fold ra n^e in angular resolution. In 'Roca-Sogorb et al.| 
(2010) we carried out a partial analysis of these data that 
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Previous 86 GHz observations of 3C 120 established an 
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Fig. 4. — Total intensity (naturally weighted) VLBA image of 
3C 120 in 2007 November 3 (2007.84) at 86 GHz. Ten contours in 
equally spaced logarithmic intervals are plotted between 0.5 and 
90% of the peak brightness of 0.54 Jy beam"-*^. The convolving 
beam of FWHM 0.37x0.27 mas at 25° is shown in the lower left 
corner. No polarization was detected at this frequency. 

revealed a brightness temperature for the component lo- 
cated 80 mas from the core a bout 600 times higher than 



expected at such distances. Roca-Sogorb et al. (2010) 
concluded that, although a helical shocked jet model may 
explain this unusually high brightness temperature, it 
appears that some other intrinsic process capable of pro- 
ducing a local burst in particle and/or magnetic energy 
may be needed to explain also the sudden appearance 
and apparent stationarity of this component. 

The 1.7 GHz total intensity image of Fig. [2] shows a 
continuous jet that extends approximately 300 mas, with 
a bulge of emission between 200 and 280 mas. Although 
the complex structure of the bulge prevents us from ob- 
taining an accurate measure of its positio n, comparison 
by eye between our image and those by | Walker et al.| 
(2001 ) - covering in total a time span of almost 30 years 
- confirms the proper mot ion of 1 . 7±0.4 mas yr~^ for the 
bulge reported by Walker et al.| (|2001|). This is similar 
to the proper motio ns measured at parsec scales (e.g., 
Gomez et al.pOOTb ), and indicates that relativistic bulk 
velocities are sustained in the jet over at least hundreds 
of par sees from the core. 

The progressively higher frequency images of Figs. [2] 
and [3] reveal a rich structure in both total and linearly 
polarized intensity even at the shortest wavelengths, be- 
ing the emi ssion at 80 mas particular ly remarkable, as 
reported by Roca-Sogorb et al. (2010). The component 
12 mas froni the core (her e after C12) is also especially 
bright. 'Roca-Sogorb et al. (2010) estimated a compres- 
sion factoij^ T] ^ 0.35 assuming that the feature corre- 
sponds to a moving shock wave. This allows the detec- 
tion of this component even at 86 GHz (Fig.B, at which 
some emission from the innermost 2 mas is also visible in 
the image. No polarization is detected at 86 GHz though. 
Comparison of the core flux density among frequencies 
yields a flat spectrum at 43 GHz and frequencies below, 
but optically thin at 86 GHz; we should note, however, 
that this may be affected by the relatively large uncer- 
tainties involved in the flux calibration at this frequency. 

^ Defined as the ratio of t he upstream to downstream particle 
densities ( Hughes et al. [19891 . 



uppe r limit to the overall core size of 54 /ias ( Gomez et al. 
1999), followed by a measurement of the FWHM size of 
the core along the axial direction of 22 ± 2 //as, as well 
as an estimation of the number density of t he combined 
electr on-positron population in the core ( |Marscher et al. 
^007|. 



3. FARADAY ROTATION ANALYSIS 

3.1. Rotation measure images from 1999, 2001, and 

2007 

Th first indications of the presence of the localized 
Faraday rotation region between 3 and 4 mas from the 
core were obtained through monthly polarimetric VLBA 
observations at 22 and 43 GHz carried out between 
November 1997 and March 1999 (Go mez et al.,,20^ 
The EVPAs of several components were ob- 



served to rotate while they passed through this region, by 
an amount which followed a quadratic dependence with 
observing wavelength - as expected for Faraday rotation. 
This is more clearly seen in the rotation measure map for 
epoch 10 January 1999 shown in Fig. [5j obtained from 
observations at 22 and 43 GHz. To solve for the ±7r am- 
biguities in the EVPAs when obtaining the RM image we 
have considered those values that result in the smallest 
RM values and yield the best agreement (both in RM 
and RM-corrected EVPAs) with the image obtained for 
observations in 2001 (see Fig. 2 of Gomez et al. 20Q8| , 
also plotted in Fig. [5] for comparison. The^^ODl KM 
image reveals the presence of an RM screen in 3C 120 
whose properties remain constant over one year, during 
which most of the jet shows rotation measure and RM- 
corrected EVPAs values with standard deviations smaller 
than 1000 rad m~^ and 20°, respectively. The Faraday 
screen displays a localized region of enhanced RM be- 
tween approximately 3 and 4 mas from the core, with a 
peak of ~ 6000 rad m~^, as well as gradients across and 
along the jet. 

Rotation measure images for observations in 2007 com- 
bining the high, intermediate, and low frequency inter- 
vals are shown in Figs. [6j [7| and [Sj respectively. 

3.2. Uncorrelated changes in jet field and Faraday 
screen 1999-2001 

To quantify any possible variation in the Faraday 
screen, we have subtracted the RM values of January 
1999 from those of the mean RM map for 2001, and com- 
puted the mean value of the difference in the components 
(jet areas in the mean 2001 map) labeled in Fig, [sl (fol - 



lowing the nomenclature of Gomez et al. 2000 2lJ01b). 



Component L is located at the localized region of high 
Faraday rotation, which is found to remain in the same 
jet area as in 2001, with very similar values of the rota- 
tion measure (mean difference of —660 ± 770 rad m~^). 
However, the RM-corrected EVPAs show significant ro- 
tation, with a mean difference between both epochs of 
—64 ± 9°. A similar situation is found for component H, 
with a mean variation in the RM of 1030 ± 1010 rad m~^ 
and a rotation of the RM-corrected EVPAs of -72 ± 19°. 
On the other hand, component K shows no significant 
changes in either the RM-corrected EVPAs or RM, with 
variations of —3 ± 19° and 810 =t 930 rad m~^, respec- 
tively. For component O we also find very similar RM- 
corrected EVPAs, with a variation of 8 ± 28°. However, 
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Fig. 5. — (top) Rotation measure map of 3C 120 in 10 January 
1999 from VLBA observations at 22 and 43 GHz. Total intensity 
(naturally weighted) 22 GHz contours are overlaid at 3.6, 5.9, 9.8, 
16, 27, 44, 73, 120, 198, 327 and 540 mJy beam"!, (bottom) Map 
of the mean value of the rotation measure for observations during 
2001 (see also Gomez et al. 2008), shown for comparison. Bars 
(of unit length) indicate the RlVl-corrected electric vector position 
angle. The lowest contour of the January 1999 image is also shown 
to facilitate comparison. 



the rotation measure in O is observed to have changed 
significantly, with a difference between both epochs of 
2730 ±910 rad m-^. 

Therefore, we find uncorrelated changes in the linear 
polarization of the underlying jet emission and the Fara- 
day rotation screen: while the RM remains constant in 
the outer jet - including the localized region of high RM 
- the RM-corrected EVPAs of two particular components 
(L and H) rotate by almost 90°; on the other hand, in 
the innermost 2 mas the RM changes significantly but 
without variations in the RM-corrected EVPAs. These 
uncorrelated changes suggest that the emitting jet and 
source of RM are not closely connected physically. 

Components L and H have RM-corrected EVPAs 
closely ali gned with the jet, i n agreement with obser- 
vations by Jorstad et al. (2005). This is contrast to what 
we find during 2001, where the RM-corrected EVPAs are 
observed to be approximately perpendicular to the jet di- 
rection. This suggests an aligned magnetic field , in agree- 
ment with that observed on larger scales (Walker et al. 



1987). If these components correspond to strong plane- 



perpendicular shock waves, compression of the magnetic 
field parallel to the shock front can result in a dominant 
mag netic field perpendicular to the components' motion 



(e.g.,'Marscher fc Gea r||1985[ [Hughes et al ||1989l [Gomez 
[et al.^ 1994), as observed in knots L and H. However, 
no sirmlar orientation of the EVPAs is found for compo- 
nent O - which is probably also associated with an even 
stronger shock - or in any of the other components that 
have traveled along the same jet regions as that of L and 
H during the 2001 observations (see also Jorstad et al. 
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Fig. 6. — Rotation measure image of 3C 120 between 15 and 43 
GHz. Contours show the 15 GHz total intensity image. Bars indi- 
cate the RM-corrected EVPAs. Inset panels show sample fits to a 
law of the EVPAs at some particular locations. The convolving 
beam is shown in the lower left corner. 



2005). On the other hand, if components L and H are 
not particularly strong shocks - which finds support by 
the fact that their degrees of polarization are similar to 
those observed for other components at similar jet loca- 
tions - we can consider the possibility that these compo- 
nents either have an intrinsically different magnetic field 
(mainly perpendicular to the jet axis), or that they have 
a different velocity and/or orientation with respect to the 
observer. In this latter case, the Lorentz transformation 
of the magnetic field - which amplifies the components 
transverse to the motion - would generally result in an 
observed dominant t oroidal field for an a ssum ed helical 
magnetic field (e.g., Blandford & Konigr"1979 Lyutikov 
et al. 2005 ) . The clear stratification in polarization ob- 
served between th e northern an d southern jet regions 
(see Fig. 3 of |GQmez et al.|[2QQ8 ) and the fact that both 
components L and H move along the southern side may 
support this idea of stratification across the jet of the 
magnetic field and/or kinematics. 

3.3. No detection of the localized high RM region 2007: 
an effect of sampling 

The middle and low frequency range RM maps (Figs.[7| 
and [8| reveal no indication of localized regions of hign 
rotation measure at these scales similar to that found 
at pa rsec scales inprevious observations ( j Gomez et al.] 
1 2 008 [ see also Fig. |5|, as may have been expected given 
the wealth of evidence for interaction between the jet and 
ambi ent mediu m (e.g.. Axo n et al.^ 1989; Sanchez e t al.| 
2004[|GQmezet al. 2000, 2008|. 'I'he high frequency range 



KM map of 1^ ig. |6] does not show the localized region of 
high RM that was found in previous observations. This 
can be understood by examining Fig. [9j which shows the 
superposition of the 15 GHz total intensity images at 
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Fig. 8. — Same as Fig. |6]for frequencies between 1.7 and 8 GHz. Rotation measure image of 3C 120 between 1.7 and 8 GHz. Contours 
show the 1.7 GHz total intensity image. 
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-5 -10 
Relative Right Ascension (mas) 

Fig. 9. — Total intensity VLB A images of 3C 120 for epochs 
2001.00 (blue) and 2007.85 (red) at 15 GHz. Both images have 
been convolved with the same beam of 1.14x0.54 mas at -1° (shown 
in the lower left corner) for easier comparison. 

the 2001.00 and 2007.85 epochs, reveahng a significant 
change in the jet geometry with time. The position an- 
gle of the innermost 3 mas structure has changed from 
-126° ± 2 in 2001.00 to -114° ± 2 in 2007.85; that is, 
the jet direction of ejection rotated to the north by 12° 
between these two epochs. As a consequence of this rota- 
tion the components in the 2007 jet do not travel across 
the localized region of enhanced RM. They thus fail to 
reveal it, since it is only through the motion of superlu- 
minal components that we are able to map the jet po- 
larization owing to the increase in energy density and 
magnetic field ordering that they produce. 

The rotation of the direction of ejectio n is in agree- 
ment w ith the helical structure suggested by [Gomez et al.| 
(|2001b|), and the proposed precessing / helical models oi 
(Japroni &: Abraham' ('2004) and Hardee et al. (2005). In 



particular, [Caproni & Abraham (2004) estimated a pre- 
cession period of 12.3 yr, predicting a swing to the north 
of the direction of ejection in the jet between our 2001 
and 2007 epochs, as is indeed observed. 

There is evidence for non-ballistic trajectories of 
the superluminal components for observations between 
November 1997 and March 1999 (Gome z et al.||2000 , 
[2001b) , resulting in the gentle bend toward the north ob- 
served in the 2001 epoch of Fig. 9] This contrasts with 
the rectilinear structure observed at the 2007.85 epoch. 
Component C12 is one of the brightest ever observed at 
parsec scales in 3C 120, so it is possible that it contains 
an unusually large energy density, with enough momen- 
tum to move ballistically and "break through" the bent 
jet funnel inferred in the 2001 epoch. This may have 
left behind a funnel through which upstream components 
move ballistically, resulting in the almost rectilinear jet 
observed in 2007.85. 

Similar changes in the geometry of the innermost struc- 
ture of the jet have also been found in the BL Lac object 
PKS 0735+178. In this case, they are related to changes 
in the overall activity of the source: quasi stationary or 
subluminal components moving through a twisted struc- 
ture are observed during periods of low activity, while 
superluminal components moving along a straight jet are 
found during periods of high activity ([G omez et al. 1999| 
2001a[ IGabuzda et al.|[2QQl| [Agudo et al.,,2006, ,Britzen[ 



et al.||2010| ) 



Figures [6] and [T] show that the rotation measure of C12 
is approximately —45 rad m~^ (— 66 ± 1 rad m~^ at the 
lower frequency interval of Fig. [8|, significantly smaller 
than the ~ 6000 rad m~^ found in the region between 
3 and 4 mas from the core (Gomez et al. 2008, see also 
Fig[5[), presumably produced by the interaction of the 
jet with the external medium. However, there is a sign 
reversal in the rotation measure of C12 when compared 
with the remainder of the jet. 

3.4. RM on larger scales from 2007 and evidence for 
reversals 

A progressive decrease in the RM along the jet in 
3C 120 is visible in the 2007 observations (Figs.[6_u| and 
[8|, as well as in the January 1999 image of Fig. bj The 
relatively high RM close to the core may be responsible 
for the dramatic decrease in the degree of polariza tion 
toward the core, as previously suggested by Gomez et al.| 
(2008). Similar decreases in the RM along the j et have 
been previousl y observe d in ot her sources (e.g., Taylor 
1998; Gabuzda & Chernetskii 20031 IJorstad et al 



laylor 



|0 'Sullivan fc Gabuzda 2009). This can be explained as 
the result ol the decrease in the density of the exter- 
nal ionized gas and magnetic field strength with distance 
from the core (e.g., iGuidetti et al.,|2010) . 

Component C12 departs from this monotonic decrease 
in RM with core distance, undergoing a change in the 
sign of the RM, which implies a double sign reversal when 
analyzing the RM along the remainder of the jet. We 
note, however, that our RM images - and therefore the 
existence of sign reversals - may be affected by a source 
of RM in the direction of 3C 120 that is unrelated to 
the jet, which recent observations estimate to be '^23 
rad m~^, with a 1-a error of ~10 rad m~^ ( [Taylor et al.] 



2009|). 

RM sign reversals have been previously reported a long 




sign reversal at parsec sc ales as reported h ere has been 
found previously (see also Guidetti et al. 2010). Sign re- 
versals in the RM across a jet can be readily explained 
by a sheath threaded by a helical magnetic field. To 
account for the reversals along a jet, different alterna- 
tive explanations have been proposed. One is based on 
"magnetic tower" models, in which the magnetic field 
lines are wound up in two nested helices with opposite 
poloidal polarities, with the inner helix less tightly wound 
than the outer helix. Such magnetic towers are expected 



to arise by the P ointing- Robert son effect (Contopoulos 
"& Kazanas"1998l), caused by the differential rotation of 
the accretion disc. A change in the relative contribu- 
tions to the observed RM of the inner and outer helical 



fields has been considered by Mahmud et al. (2009) to 



explain changes in the direction ol transverse KM gra- 
dients, which can eventually also lead to a sign reversal. 
Similar magnetic tower models have also been consid- 
ered to explain the different helical pitch angles derived 



for the emitt ing jet and the sheath i n NRAO 140 ( Asada 
et al.[|2008b[ ) and S5 0836+710 (As ada et al.|[2010] ). It is 



II component C12 is interacting with the external 
medium, it does not result in enhanced rotation measure. 



also possible to explain a sign reversal in the KM along 
the jet by assuming that the sheath responsible for the 
Faraday rotation is at least moderately relativistic, so 
that light aberration can lead to a change in the line-of- 



On the source of Faraday rotation in 3C 120 



9 



sight dominant component of the magnetic field through 
a small change in the dire ction or velocity of the flow, a s 
considered previously by p'Sullivan fc Gabuzd a' f2QQ9V 
It would be difficult to explain the double RM sign 
reversal along the jet of 3C 120 trough magnetic tower 
models, since this would consequently require an unlikely 
double change in the relative contributions to the RM of 
the inner and outer helical flelds. However, the change in 
the RM sign in C12, when compared with the RM of the 
remaining jet, can be explained within the moderately 
relativistic sheath model as produced by a local change 
in the velocity or orientation, as would be expected in the 
case of precession in 3C 120. Note that this could also 
explain its unusually bright emission through differential 
Doppler boosting. 

It is interesting to note that, although the RM of C12 
undergoes a sign reversal, the RM-corrected EVPAs have 
the same orientation as for the remainder of the optically 
thin jet, i.e., perpendicular to the jet. (Note that com- 
ponent C80 contains RM-c orrected EVPAs aligned with 
the jet, but as reported by |Roca-Sogorb et al.||2010[ this 
is probably due to its peculiar nature.) Therefore, in 
order to explain the RM sign reversal in C12 through 
any of the different proposed models without produc- 
ing a change in the RM-corrected EVPAs, we need to 
assume that the emitting jet and the sheath responsible 
for the Faraday screen have different kinematical proper- 
ties and/or magnetic fleld conflgurations. Alternatively, 
a much simpler interpretation for the different RM sign 
in C12 can be obtained by considering that the Fara- 
day rotation is produced in a foreground cloud with a 
different magnetic fleld than that of the remaining jet. 

A sign reversal in the rotation measure alongthe jet is 
also observed in the January 1999 image (Fig.|5|, as well 
as between the innermost 2 mas of t he jet in this image 



and those obtained in 2001 (Fig. 1 in Gomez et aL[2008 ) 
and 2007.85 (Fig. |6|. If we assume that the sheatn is at 
least moderately relativistic, a change in the direction of 
ejection in the jet - as expected for precession - could 
explain the differences in the RM between 1999 and 2007 
through the different projection of the fleld along the line 
of sight. To explain the RM sign reversal between 1999 
and 2001, an unlikely large change in the bulk velocity 
of the sheath would be needed (given the similar direc- 
tions of ejection for the jet at both epochs), so that a 
change in the line-of-sight component of the fleld could 
be obtained through light aberration. Note however that, 
as mentioned for C12, different kinematical properties 
and/or magnetic fleld conflgurations would be needed for 
the sheath and emitting jet, so that the RM can change 
without producing a rotation of the RM-corrected EV- 
PAs^ 

Broderick & Loeb" f2009) have considered the possi 



bility of having relativistic helical motions in the Fara- 
day rotating sheath, showing that sign reversal in the 
rotation measure can naturally arise in this case, espe- 
cially for observations at high frequencies that probe the 
innermost regions where there is a transition between 
ultra-relativistic and moderately relativistic helical mo- 
tion. This could explain the RM sign reversals along the 
jet in the January 1999 image, but it would require a - a 
priori unlikely - considerable change in the kinematical 
properties of the sheath to account for the RM time vari- 
ability observed between January 1999 and later epochs. 



The relat ivistic helical sheath model of lBroderick fcl 
Loeb (2009) might also account for the sign reversal of 
component C12, which could also explain the variation 
in the RM with changing frequency intervals. 

4. SUMMARY AND CONCLUSIONS 

We have studied the source of the Faraday rotation 
screen in the radio galaxy 3C 120 through VLB A obser- 
vations carried out between January 1999 and November 
2007. Sampling of the rotation measure from the inner- 
most regions up to hundreds of mas from the core is 
obtained in November 2007 through observations at 86, 
43, 22, 15, 12, 8, 5, 2, and 1.7 GHz, providing a 50-fold 
range in angular resolution. 

By analyzing correlated motions of some particular su- 
perluminal components when performing the registration 
of images across epochs, we show that it is possible to 
achieve an accuracy between 12 and 30 /ias, similar to 
what can be obtained by phase-reference VLBI astrom- 
etry. 

Signiflcant changes are observed in the RM when com- 
paring observations taken in 1999 and 2001 with those 
in 2007. In particular, one of the main features observed 
during the 1999 to 2001 observations, the existence of a 
localized, very stable high RM region between 2 and 3 
mas from the core, is not observed. This is explained not 
as a result of an actual change in the Faraday screen, but 
by a transverse shift in the positions of bright features in 
the jet, such that the screen no longer covers superlumi- 
nal knots 2-3 mas from the core. A similar explanation 
has been off ered to explain the RM variability in 3C 279 
and 3C 273 ( [Taylor |20Q0| [Zavala fc Taylor |200lD . The dif- 
ferent trajectories of the super luminal components (i.e., 
different jet geometry) are in agreemen t with p roposed 
precession models f or 3C 120 (Caproni fc Abraham 2004^ 
Hardee et alT][2QQ5[ ). 

Comparison of observations taken between January 
1999 and December 2001 reveals uncorrelated changes in 
the linear polarization of the underlying jet emission and 
the Faraday rotation screen: while the rotation measure 
remains constant between approximately 2 and 5 mas 
from the core during these three years of observations, 
the RM-corrected EVPAs of two particular superluminal 
components (moving along the southern side of the jet) 
are found to be rotated by almost 90° when compared 
to other components moving through similar locations 
in the jet. On the other hand, in the innermost 2 mas 
the RM changes signiflcantly - including a sign reversal 
- but shows no variations in the RM-corrected EVPAs. 
These uncorrelated changes suggest that the emitting jet 
and the source of the RM are not closely connected phys- 
ically. 

Similarly, observations in 2007 reveal a double sign re- 
versal in the RM along the jet while the RM-corrected 
EVPAs remain perpendicular to the jet axis. The varia- 
tion in the RM can be accounted for by considering that 
the sheath responsible for the Faraday rotation is at least 
moderately relativistic, so that light aberration can lead 
to a reversal in the polarity of the line-of-sight magnetic 
fleld through a small change in the direction and/or ve- 
locity in the flow of the sheath. Note, however, that dif- 
ferent kinematical properties and/or magnetic fleld con- 
flgurations would be needed for the sheath and emitting 
jet, so that the RM can change without producing a ro- 
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tation of the RM-corrected EVPAs. Such differences in 
magnetic field structure can be explained by a magnetic 
tower model, in which the inner helix would thread the 
emitting jet and the outer helix the RM sheath, but an 
unlikely double change in the ratio of the RM of the in- 
ner to the outer helix would be required to explain the 
double sign reversal in RM along the jet. 

The observed coherent structure and gradient of the 
RM along the jet supports the idea that the Faraday ro- 
tation is produced by a sheath of thermal electrons that 
surrounds the emitting jet. However, the uncorrelated 
changes in the RM screen and RM-corrected EVPAs sug- 
gest that the emitting jet and the source of Faraday rota- 
tion are not closely connected physically and have differ- 
ent configurations for the magnetic field and/or kinemat- 
ical properties. Furthermore, the existence of a three- 
year-long stationary region of enhanced RM requires a lo- 
calized source of Faraday rotation, which favors a model 
in which a significant fraction of the Faraday rotation 
measure found in 3C 120 originates in foreground clouds, 
rather than in a sheath intimately associated with the 
emitting jet. In this case, Faraday rotation studies will 
provide valuable information about the ambient medium 



through which jets propagate, but will not be able to re- 
veal further details about the emitting jet, such as the 
line-of-sight magnetic field, and hence to test whether 
they are threaded by helical magnetic fields. 
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APPENDIX 



REGISTERING OF THE 2001 IMAGES 

Analysis of the images corresponding to the 2001 monitoring (Fig. [T]) across epochs and frequencies requires first 
a proper alignment (or registration) of the maps, given that self-calibration of the visibility data during the imaging 
process results in the loss of absolute positioning of the source. To avoid this it is possible to perform phase-referenced 
observations, in which the position of the target source can be determined relative to a reference source. However, this 
requires the reference source to be strong, close to the target source, and ideally point-like. These conditions are rarely 
met, and therefore the alignment of VLBI images is commonly performed by looking for particular jet feature(s) that 
can be unambiguously identified in each of the images to be compared. The core (defined here as the position where the 
jet becomes optically thin) seems a priori a good choice, since it is easy to identify. However, opac ity effects can move 
its position downstream as the frequency of observation decreases (e.g., iBlandford & Konigl 1979 Gomez, Alberdi, & 
Marcaide|1993 Lobanov|1998 Caproni fc Abraham|2004 Hirotani|2005 [Kovalev et a l. 2008), rendering the alignment 
of images at multiple frequencies a ditticult task. I^'urtherniore, the ejection of a new component produces a temporary 
dragging of the apparent core position downstream. This is due to the increased opacity and the motion of the peak 
brightness resulting from the convolution of the actual core and the new moving component. Variations in the jet 
flow velocity or vi ewing angle, as expected for precessing jets, should also lead to a change in the core position (e.g., 
Caproni fc Abraham 2004). This jitter of the core position complicates the alignment of images, especially in active 
sources like 3 (J 120. It is therefore more convenient to choose components that are known to be optically thin, so that 
opacity effects do not affect the position when images are compared at different frequencies. To minimize the errors 
in the determination of the position such components are ideally strong and compact, that is, have high brightness 
temperatures. Finding such components is not always easy, since jet features can be observed to split or merge , or 
to be diff i cult t o disentangle from other neighboring components. In an effort to overcome such difficulties, Cro ke &: 
Gabuzda (2008) describe a new method based on ima ge cross-correlation tec hniques, and show that it can be used 
to align maps at different frequencies. More recently, jPodson fc Rioja (2009) have reported a novel phase-reference 
procedure, dubbed source /frequency phase referencing (Sl^PK), which combines the source switching of conventional 
phase-referencing techniques with fast frequency switching, allowing astrometric measurements of frequency-dependent 
core positions up to 86 GHz. 

The characterization of the images of Fig.[l]in terms of discrete jet components has been carried out using DIFMAP, 
which allows one to fit the visibility data with Gaussian brightness distributions. Multi ple jet components can be 



identified in the images of Fig. [Tl Most of these display complex e volution (see also, e.g., Gomez et al. 2000 2001b 



Homan et aT]|2001[ [Jorstad et aL]|2005[ |Agudo et aL]|2007| [2010'), or are difficult to match at the three observms; 



frequencies simultaneously, and therefore cannot be used to register the images. However, components 6, c, and / at 
15 GHz, and b and / at 22 and 43 GHz (see Fig.[l) can be very reliably identified across epochs and frequencies, and 
therefore can be used to align the maps. We refernereafter to these as the "registering components" . 

DIFMAP sets the coordinate origin at the phase center of the visibility data, that is, the peak in the initial map. 
This is always close to, but not necessarily coincident with, the peak brightness of the final clean image, which in 
turn usually coincides with the core as long as the jet is core dominated. That is, the origin of coordinates is usually 
coincident with the core. But as we have discussed above, the core position depends on the frequency of observation 
and can move when new components are being ejected. Therefore, the positions of the components are actually referred 



On the source of Faraday rotation in 3C 120 



11 






2000.8 2001.0 2001.2 2001.4 2001.6 2001.8 2002.0 
Epoch (yr) 



2000.8 2001.0 2001.2 2001.4 2001.6 2001.8 2002.0 
Epoch (yr) 



2000.8 2001.0 2001.2 2001.4 2001.6 2001.8 2002.0 
Epoch (yr) 




2000.8 2001.0 2001.2 2001.4 2001.6 2001.8 2002.0 
Epoch (yr) 




1.0 2001.2 2001.4 2001.6 2001.8 2002.0 
Epoch (yr) 




2000.8 2001 .0 2001 



2 2001 .4 2001 .6 2001 
Epoch (yr) 



Fig. 10. — Residuals of the linear fits in right ascension and declination with time for the registering components at 15 GHz (left), 22 
(middle) and 43 GHz (right). Blue triangles correspond to component b, green asterisks to c and red crosses to /. Component c appears 
blended with other jet features at 43 GHz in epoch A, and therefore has not been used in the analysis. 

to a location that is subject to change with epoch and frequency, and by an amount that it is unknown. However, 
there is a possible way to at least estimate these core motions. All jet components are referred to the core position, 
and therefore any motion of the core through time or frequency should affect all of them by the same amount. To 
search for these correlated motions we have considered the simplest case of constant rectilinear motion, which in fact is 
an excellent representation of the actual motion of the registering components, at least during the time covered by our 
2001 observations. Core motions should appear as correlated deviations of the components from these simple linear 
fits. We should note, however, that linear motions of the core itself would obviously not be revealed by this procedure. 

Figure [lO] plots the residuals of linear fits in right ascension and declination for the registering components at the 
three frequencies. We have excluded the data for the last epoch in this analysis because of their poor quality. Figure 10 
clearly shows a pattern, far from the random behavior that would be expected in the case of changes in position caused 
by systematic errors in the model fitting, or by changes in the internal structure of the components. Accelerations, 
including bent trajectories, would also produce significant residuals, with a clear time evolution, as shown in Fig. \iO\ 
However, it is highly unlikely that such deviations would be the same and occur simultaneously at two (three in 
the case of 15 GHz) components that are located in very different jet regions. (Components b and c move between 
approximately 1 and 3 mas from the core, while component / moves between about 6.5 and 8 mas.) We therefore 
conclude that the residuals of Fig. 10 are representative of the core motion across epochs. To quantify this motion 



we have used the mean value of the residuals at each epoch, with an estimated error given by the standard deviation. 
The mean values of the errors are 15, 14, and 12 /ias in right ascension, and 22, 21, and 30 /ias in declination, at 15, 
22 and 43 GHz, respectively. 

This procedure therefore provides an acc uracy in the a lignment of the im ages that is similar to w hat can be achieved 
by phase-reference VLBI astrometry (e.g., iGuirado et al. 2000) or SFPR (Dodson & Rioja 2009). Alignment of the 
images across frequencies is obtained through matching of the position of the optically thin registering components c 
and /. 
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